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1. Introduction 

Staggered fermions are computationally faster than Wilson fermions, the discretization errors 
scale with 0'{a^) and due to the well behaving spectrum of the Dirac operator, light quark masses - 
nowadays even the physical ones - can be reached. Note, however, that a straightforward definition 
for staggered fermions exists only for 4,8,... quark flavours. All works using 2 or 2 + 1 flavours of 
staggered quarks use the fourth root (or square root) trick to have only one (or two) flavour(s). The 
action is defined by taking the fourth root of the fermion determinant. Since such a prescription is 
non-local at fixed lattice spacings, it is debated whether it is equivalent to a local continuum field 
theory or not (see e.g. [1-3] and references therein). 

Large scale staggered studies are based on the expectation that in the continuum limit staggered 
and Wilson results agree. Clearly, at fixed, non-vanishing lattice spacings deviations are possible 
even if the continuum results are the same. Therefore, any analysis looking for the equivalence or 
non-equivalence between the staggered and Wilson formalism should fulfill two conditions. First 
of all, the analysis should be based on several lattice spacings and a controlled continuum ex- 
trapolation should be carried out. At least three different lattice spacings are needed, for which the 
asymptotic scaling behaviour can be already observed (both for staggered and for nonperturbatively 
clover-improved Wilson fermions one expects an a^ scaling). Secondly, the physical quantity cho- 
sen as a basis of such comparison not only has to be well defined and relatively easily measurable, 
but also has to be sensitive to the dynamical fermion sector. 

It is not so obvious how to find such a sensitive quantity. E.g. the masses of different hadrons 
provide well defined physical quantities, however, these are not sensitive enough to the dynamical 
fermion sector. Even the quenched calculations, where the dynamical fermions are completely 
omitted, provide the physical hadron masses with an error about or less then 10% [4]. Finding 
measurable differences in the continuum limit of such quantities would require extremely high 
precision, thus, prohibitively large scale calculations. 

Thermodynamic observables can be much more sensitive to the fermionic content of the the- 
ory. This fact is related to the singular/non-singular behaviour of the finite temperature QCD tran- 
sition. E.g. on the one hand the quenched theory undergoes [5, 6] a first order phase transition at 
non-vanishing temperatures (T). For first order phase transitions different observables behave in a 
singular way (the latent heat is infinite or the temperature derivative of the renormalized Polyakov 
loop has a discontinuity). On the other hand the finite temperature transition of QCD with staggered 
fermions and physical quark masses has turned out to be a crossover [7]. As a consequence, none 
of the physical quantities as the function of the temperature are infinite or discontinuous, at the 
most they only undergo a rapid change within a narrow temperature range. Clearly, the height of 
such a peak is expected to be quite sensitive to the details of the fermionic properties of the action. 
As we increase or decrease the quark masses the peak turns out to be more and more singular and 
after a while a second order then a first order phase transition region is reached. Thus, the fermionic 
content (quenched or unquenched with physical or non-physical quark masses) manifests itself in 
a very pronounced way. 

Choosing a physical quantity that undergoes a rapid change and finding the maximum of its 
derivative with respect to the temperature may provide a quantity sensitive to the finite temperature 
behaviour of the system. 
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The aim of this study is to provide a quantitative comparison between the behaviour of the 
staggered and Wilson fermion formulations. We attempt to perform an analysis at two different 
lattice spacings (thus no conclusive continuum extrapolation can be done yet). The observable we 
have chosen is related to the transition at non-vanishing temperatures and expected to be sensitive 
to the fermionic properties of the theory. 

One physical quantity that undergoes a rapid change around the transition temperature is the 
quark number susceptibility, which is defined via [8] 



Zv 1 5^1ogZ 



r2 TV d^i^ 



(1-1) 

;U=0 



where /x is the quark chemical potential. The quark number susceptibility can be directly mea- 
sured, and it automatically has the correct continuum limit without the need for renormalization. 
In addition, the maximum of its derivative with respect to the temperature, that is, the rate at which 
the susceptibility changes during the transition, is sensitive to the dynamical fermion sector. These 
properties make the quark number susceptibility a good candidate for the quantity that a comparison 
of the Wilson and staggered fermion formulations should be based on. Since the rooting procedure 
of the determinant is less transparent (might be more problematic) for odd number of flavours we 
use three flavours. (The one flavour theory, where there is no chiral symmetry breaking [9], is not 
suitable for our purposes.) 

2. Action parameters 

The three flavours were degenerate for both Wilson and staggered calculations, and the gauge 
action used was the Symanzik tree-level improved gauge action [10]. In both cases the lattices 
sizes were 2)2^ x 8 and 32^ x 10, and the configurations were generated using the Rational Hybrid 
Monte-Carlo algorithm [11]. 

Wilson calculations: Three steps of stout smearing [12] with smearing parameter p = 0.1 were 
used. The gauge coupling constant was in the range j3 = 3.2 — 3.7. In addition the femionic sector 
was clover improved [13] with a tree level clover coefficient c = 1.0. Note, that for this type of 
smeared fermions the tree level clover coefficient essentially leads to an ^(a) improved action [14]. 

Staggered calculations: Two steps of stout smearing with smearing parameter p =0.15 was 
used [15], and the gauge coupling constant was in the range j3 = 3.5—4.0. This staggered action 
within this lattice spacing range was shown to be in the scaling regime [16]. 

3. Setting the scale 

To be able to make sure that the staggered and Wilson calculations are performed at the same 
set of physical parameters, the line of constant physics was defined by fixing the ratio of the pseu- 
doscalar and the vector meson masses mps/my. 

If the finite temperature behaviour of the Wilson and the staggered formulations are different, 
the difference is likely to be most apparent at small quark mass parameters. For small quark masses 
we are closer and closer to the first order phase transition region, thus all the differences are easier 
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Figure 1: The scale for the staggered calculations: the bare staggered quark mass (left panel) and the vector 
meson mass (right panel) as a function of the gauge coupling constant. 
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Figure 2: The scale for the Wilson calculations: the bare Wilson quark mass (left panel), the PCAC quark 
mass (middle panel) and the vector meson mass (right panel) as a function of the gauge coupling constant. 



to see in the height of the peak of our observable. Decreasing the quark mass, however, causes 
the computational costs to rise. Therefore, one needs to find a compromise between making the 
comparison more sensitive by lowering the quark mass (or in other words lowering nips/mv) and 
keeping the computational costs reasonable. 

As a compromise we have chosen mps/my = 0.5, which sets the quark mass about m,,./3, 
where my is the physical strange quark mass. This relationship completely defines the line of 
constant physics. In three flavour lattice QCD we have two parameters. One of them is the quark 
mass which is essentially set by the relationship between the pseudoscalar and vector mass ratio. 
The other one is the lattice spacing, which is dominantly given by the gauge coupling. 

The bare quark mass and the vector meson mass corresponding to the different gauge couplig 
values for the staggered case are shown in Figure |l]. For the Wilson calculations, the bare quark 
mass, the current algebra quark mass and the vector meson mass are shown in Figure 0. 
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Figure 3: The quark number susceptibility (left panel) and its derivative (right panel) as a function of the 
temperature normalized with the vector meson mass. 

4. Results 

To be able to compare the staggered and Wilson results the temperature was made dimension- 
less by dividing by the vector meson mass. The quark number susceptibilities as a function of 
T jmy are shown in the left panel of Figure ^. The derivative of the susceptibility was obtained by 
fitting cubic polynomials to the susceptibility points, then taking the derivative of the polynomial. 
The slight change due to the variation of the fitting range is taken as a systematic error. These 
derivatives are shown in the right panel of Figure ^. 

5. Conclusions 

Calculating the maximum of the derivative of the quark number susceptibility with respect to 
the temperature yields a sensitive quantity, based on which the finite temperature behaviour of the 
Wilson and the staggered fermion formulations can be compared. We have performed calculations 
using lattices with temporal extensions Nt = S and 10. At A'f = 8 the maximum of the dervative in 
the Wilson case is a factor of 2 higher than the staggered result, whereas at A'; = 10 the heights of 
the peaks get closer to one another. In order to be conclusive one needs to get one or more steps 
closer to the continuum limit by performing calculations using Nt = 12 or even finer lattices. 
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